Time Delay in Black Hole Gravitational Lensing as a Distance Estimator 
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We calculate the time delay between different relativistic images formed by black hole gravitational 
lensing in the strong field limit. For spherically symmetric black holes, it turns out that the time 
delay between the first two images is proportional to the minimum impact angle. Their ratio gives 
a very interesting and precise measure of the distance of the black hole. Moreover, using also the 
separation between the images and their luminosity ratio, it is possible to extract the mass of the 
black hole. The time delay for the black hole at the center of our Galaxy is just few minutes, but 
for supermassive black holes with M = 10 s -f- 1O 9 M0 in the neighbourhood of the Local Group the 
time delay amounts to few days, thus being measurable with a good accuracy. 
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I. INTRODUCTION 

Gravitational lensing is a useful tool to investigate a lot 
of aspects of the nature of the universe. It was the first 
prove of the validity of the theory of general relativity 
(GR) 0, and today its effects on extragalactic scales 
(lensing of quasars, arcs in galaxy clusters, etc.) and 
on galactic scales (microlensing) are ordinarily observed 
and studied by the scientific community in the weak field 
approximation [2|. 

In the last years, a new form of gravitational lensing 
has been proposed as a method to investigate the gravita- 
tional field generated by collapsed objects. This approach 
considers light rays of background sources passing very 
close to the event horizons of black holes without entering 
inside. The study of this extreme case is of remarkable 
interest: on one hand it represents an independent test of 
GR in strong gravitational fields; on the other hand, in- 
trinsic features of the lens (rotation, electric charge, etc.) 
could become accessible to the observations, opening a 
new possibility to constrain black hole models. 

In Schwarzschild framework, a light ray with small im- 
pact parameter can wind several times around a black 
hole without being definitively caught inside. In this 
way, a set of infinite relativistic images can be gener- 
ated on each side of the black hole y, 0- Several ap- 
proaches have been proposed to study gravitational lens- 
ing in the strong field limit 0, H ll S HH E3 , while 
Falcke, Melia & Agol, in a different perspective, stud- 
ied the accretion flow as a source [llj. In Ref. 01 
analytical formulae for the position and the magnifica- 
tion of the images were obtained, defining a strong field 
limit for the deflection angle. These formulae were ap- 
plied to a Reissner-Nordstrom black hole in Ref. |13| . 
and were also used to calculate relativistic effects on mi- 
crolensing events 0] . A full generalization of the strong 
field limit for any spherically symmetric spacetime was 
drawn in Ref. 15] and applied to several black hole met- 



rics, allowing a non-degenerate discrimination among dif- 
ferent collapsed objects. In work 0, the method was 
used to examine the characteristics of a Gibbons-Maeda- 
Garfinkle-Horowitz-Strominger charged black hole of het- 
crotic string theory. In work 0] the case of a Kerr black 
hole and the relevance of its spin in strong field lensing 
approximation was discussed for light rays travelling on 
quasi-equatorial trajectories. Waiting for an analytical 
treatment including non-equatorial trajectories, the gen- 
eral case is explored numerically in |18| . 

When multiple images are formed, the light-travel-time 
along light paths corresponding to different images is gen- 
erally not the same. So, if the source is characterized 
by luminosity variations, these variations would show up 
in the images with a relative temporal phase depending 
on the geometry of the lens [l!j . These time delays are 
usually measured in gravitational lensing observations on 
cosmological scales. The striking importance of time de- 
lay lies in the fact that it is the only dimensional observ- 
able. Therefore its measurement is useful to determine 
at first the length scale for a gravitational lensing system 
and its mass. Measuring the time delays in cosmological 
contexts, it is possible to determine the cosm olog ical dis- 
tance scale and hence the Hubble parameter |l9ll20t l2l| . 
This fact has drawn a great attention by the scientific 
community towards this kind of measurements. 

In the present paper, we estimate the time delay be- 
tween images generated by strong field lensing of black 
holes. We show that time delays between relativistic im- 
ages are indeed measurable in most supermassive black 
holes suitable for gravitational lensing studies. Moreover, 
it turns out that in a first approximation the time delay 
between consecutive relativistic images is proportional to 
the minimum impact angle. The ratio between these two 
observables is nothing but the distance of the lens, that 
can be estimated in a very precise way and without bias. 
Combining all information, it is also possible to get an in- 
dependent mass estimate, to be compared with estimates 
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obtained by other methods. 

This paper is structured as follows. In Sect. 2, we re- 
call the main results of the strong field limit method. In 
Sect. 3 we derive a general expression for the time delay, 
specifying it to the spherically symmetric case. In Sect 
4, we estimate the expected time delays for several inter- 
esting supermassive extragalactic black holes discussing 
the whole information that can be extracted from a time 
delay measurement. Finally, in Sect. 6, we draw the 
conclusions. An appendix contains the computation of 
the time delay for the Kerr metric, as an example of a 
non-spherically symmetric metric. 



II. THE STRONG FIELD LIMIT APPROACH 

The technique we use in the derivation of the time 
delay between different images resembles the main cal- 
culation of the deflection angle in the strong field limit 
approach. We shall briefly recall the main steps of that 
derivation referring the reader to Refs. for all the 

details. 

Consider a generic black hole metric projected on the 
equatorial plane 

ds 2 = A(x)dt 2 - B(x)dx 2 - C{x)d<t> 2 + D{x)dtd(f> (1) 

where D(x) can be consistently set to zero in a spherically 
symmetric black hole. 

The metric does not depend on time and the azimuthal 
angle (f>, so that, for a photon moving in this background, 
t and 4> can be expressed in terms of two integrals of 
motion, namely energy and angular momentum. By a 
suitable choice of the affine parameter, we set the first 
to 1 and the second to the impact parameter u of the 
incoming photon. We have 



4C - 2uD 

4AC + D 2 
4Au + 2D 

= 4AC + D 2 ' 



(2) 
(3) 



The impact parameter u, is related to the closest ap- 
proach distance Xq by 



-Do + y/4A C + D 2 

2An 



(4) 



where all functions with the subscript are evaluated for 

x = Xq. 

By the on-shell condition for the photon, we also derive 
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VBV 4AC + D 2 
Dividing Eq. J3J by Eq. JSJ), we get 



(5) 



d4 



P 1 (x,x )P 2 (x,x Q ) 



Pi(x,x ) 



P2(X,X ) 



B{2A Au + A a D) 

^/U^VaacTd 2 
i 



A - A%- + %(AD - A D) 



(6) 
(7) 
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Integrating this expression from xq to infinity we find 
half the deflection angle as a function of the closest ap- 
proach. Given the symmetry between approach and de- 
parture, we can write the whole deflection angle as 

(9) 
(10) 



a(x ) = <pf(x ) - 7T 
M*o) = 2j -dx. 



To solve this integral, we define the variables 

y = A{x) (11) 

z = ^ (12) 
l-2/o 

where yo = A)- The integral llUfl in the deflection angle 
becomes 

l 

^>f{xo) = J R{z,x )f(z,x )dz (13) 
o 



R(z,x ) = 2 1 V " Pi{x,x ) 
A'{x) 



(14) 

f(z,x ) = P 2 {x,xo) (15) 

where x = A^ 1 [(1 - y Q ) z + y }. 

The function R{z, x$) is regular for all values of z and 
xq, while f(z,xo) diverges for z — > 0. We then expand 
the argument of the square root in f(z, Xq) to the second 
order in z, defining 

f(z, x ) - f (z, x ) = , 1 n „ ■ (16) 

az + pz l 

The Eq. a = defines the radius of the photon sphere 
x m , which is the minimum approach distance for photons 
not falling into the black hole. 

The result of the integral (|13|) gives the strong field 
limit expansion of the deflection angle 



(17) 

(18) 
(19) 
(20) 



a{u) = —a log 1 +b + O (u — u m ) 

where the coefficients of the expansion are 



D m + ^AA m C m + D 2 m 



2A n 



R{0,x m ) 

a = ■ 

2^ 



it + b D + b R + alo; 
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and 



b D = 2a log 



2(1 - y m ) 

A' m X m 



(21) 



bn = J [R{z,x m )f(z,x m ) - R(0,x m )fo(z,x m )]dz, (22) 
o 

while c is defined by the expansion 

u - u m =c(x - x m ) 2 . (23) 
All the functions with the subscript m are evaluated at 

Xq X nl . 

With the formula l|17fl for the deflection angle, it is 
straightforward to calculate the positions and the mag- 
nifications of all relativistic images. Two infinite patterns 
of relativistic images appear on each side of the lens, very 
close to the minimum impact angle 9 rn — u m /DoL (Dol 
is the distance of the lens from the observer). These 
images are highly demagnified unless the source is very 
close to a caustic point. For spherically symmetric black 
holes, all caustic points are exactly aligned with the lens, 
so that a source aligned with the optical axis (the line 
joining observer and lens) would enhance the magnifica- 
tion of all images simultaneously. 

In spinning black holes, the caustics drift away from 
the optical axis, so that one source cannot be simulta- 
neously close to different caustics. In this case only one 
image at a time can be enhanced while all others stay 
very faint 01- Nevertheless, in this case, additional im- 
ages, appearing when the source is inside a caustic, may 
play an important role in the phenomenology, yet to be 
understood. 

For later reference, we write here the formula for the 
position of the relativistic images 



= ±9 m 1 



T±- 



(24) 

Here 7 is the angular separation between the source and 
the optical axis, as seen from the lens, n is the number of 
loops done by the photon around the black hole. For each 
n, we have an image on each side of the lens, according 
to the chosen sign. 



III. TIME DELAY IN THE STRONG FIELD 
LIMIT 

In this section we derive the time delay between differ- 
ent relativistic images, following an approach similar to 
the one reported in the previous subsection for the deflec- 
tion angle, but with some tricky subtraction strategies to 
treat the integrals. 

For an observer at infinity, the time taken from the 
photon to travel from the source to the observer is simply 



dt. 



(25) 



Changing the integration variable from t to x, we split 
the integral into approach and leaving phases 



xo Dol 

T = / —dx + / —dx. 

J dx J dx 



(26) 



Here Dls is the distance between the source and the 
lens, while Dol is the distance between the lens and the 
observer. 

Extending the integration limits to infinity, we can 
unify the two integrals into one, exploiting the symmetry 
between approach and departure. This can be done at 
the price of subtracting two terms 



T = 2 



dt 
dx 



dx — 



dt 
dx 



dx 



dt 

dx 



dx. (27) 



If we consider two photons, travelling on different tra- 
jectories, the time delay between them is 



Ti -T 2 = 2 



dt 



Dol 



D l 



dt 
dx 



(x,x 0A ) 



dt 

Tx {w) 



dx 



dx 



dx-2 



dt 



dt 



dx 



Dol 



D l 



dx 



dt 



dx 



dx. (28) 



Supposing that observer and source are very far from 
the black hole, dt/dx is effectively 1 in the last four inte- 
grals which thus exactly cancel each other. We are thus 
left with the first two integrals. 

Dividing Eq. J2J by Eq. (JjJJ, we obtain 



= Pi(x,x a )P 2 (x,xo) 

dx 



Pi{x,Xq) = 



__ y/BM{2C - uD) 



(29) 
(30) 



to 



and P2 defined by Eq. (JHJ. Of course, dt/dx tends to 
one for large x and the two integrals in (|28|l are sepa- 
rately divergent, while their difference is finite. In fact, 
the time delay is the result of the different paths followed 
by the photons while they wind around the black hole. 
When the two photons are far away from the black hole, 
dt/dx — > 1 and the two integrals compensate each other. 
Separating the two regimes, we can write individually 
convergent integrals. To achieve this, we subtract and 
add the function Pi(x, Xq^)/ -J Aq^ to each integrand. 
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Supposing xq.i < X0 2, we can write 



T 1 -T 2 = T(x 0A ) - T(x 0}2 ) + 2 



x 0,2 _ 

Pi(x,x ,i) 



dx 



m n S 
Q- 



with 



Pi (21,3:0,1) Pi(aCj»o,2) 



fi;j; 



(31) 



S 

--Q- 



f(x )= / R(z, x )f(z,x )dz 



R(z,x ) 



A'(x) 



Fi{x,X ) 



1 - 



1 



fAof(z,x ) 



(32) 



(33) 



and f(x,Xo) defined by Eq. I|15[). Substituting all the 
expressions back into (|31|l . we can check that it is equiv- 
alent to (|28|l . but now it is written as a sum of separately 
convergent integrals. 

In practice, the integral T(xq) represents the time 
spent by the light ray to wind around the black hole. 
In order to cutoff the integrands at large x's, in the def- 
inition of R(z,xq) we have subtracted a term which is 
negligible when the photon is close to the black hole but 
cancels the integrand when the photon is far from the 
black hole. The residual terms of this subtraction are 
stored in the last two integrals in (|31|l and are generally 
subleading with respect to AT, as we shall see later. 

The integral (|32|l can be solved following the same tech- 
nique of the integral 113|) in the previous subsection, just 
replacing R by R. The result is 



(34) 



T(u) = -5 log 1 + b + O (u- u m ) 

where u m is defined by Eq. 118f) and 

. R(0,X m ) 
a = , 

2v^ 



with 



bo = 2a log 



b = -7r + b D + b R + a log ■ 



2(1 - y m ) 
A' m x m 



(35) 
(36) 

(37) 



b R = 



R{z, x m )f{z, x m ) - R(0, x m )f {z, x m ) dz. (38) 



For spherically symmetric black holes, the expression 
for the time delay can be advantageously simplified. No- 
tice that, for spherically symmetric spacetimes, D = 
and 



Pi(x,xo)|_d = o 



BAn 



(39) 



n S 
G- 



FIG. 1: This figure illustrates which images are considered in 
the different time delay calculations in the text. L represents 
the lens, S is the source. On each side of the lens an infinite 
series of images is formed. In case a, we consider two images 
on the same side of the source. Their time delay is given by 
AT^ im with the upper sign. In case 6, we consider two images 
appearing on the same side but opposite to the source. Their 
time delay is given by AT° m with the lower sign. Finally, in 
case c, we consider two images appearing on opposite sides. 
Their time delay is given by AT° m . 



Then the last integral in (|31|l identically vanishes. When 
D 7^ 0, the dependence on xq remains through the impact 
parameter u which is present in Eq. (|30|l . The second 
integral in l|3I(l can be approximated substituting the in- 
tegrand with ^/B m /A m since it is practically constant 
throughout the (very small) integration interval. Finally, 
combining (|34() with (|17f) we get a very simple expression 
for the first term in l|31|) . 

In writing the final formulae, we distinguish the case 
when the two images are on the same side of the lens 
from the case when the two images are on opposite sides 
of the lens. In the first case, we have 



AT* = 2tt(u 




(40) 



where the upper sign before 7 applies if both images are 
on the same side of the source (Fig. [TJi) and the lower 
sign if both images are on the other side (Fig. ^>). 

If the images are on opposite sides of the lens (Fig. 
lit), then 



ATZ m = IMn - m) - 2 7 ] - 



„ / B m ju m b_ ( 



,(41) 



where the image winding n times is on the same side of 
the source and the other is on the opposite side. 
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Notice that the geometry mostly favoured for the ob- 
servation of relativistic images is that with the source 
almost aligned with the lens, so that 7 
Therefore 



D OL « 2tt. 



n,n ~~ n,m^n 



rpO 



(42) 



i.e. if we evaluate the time delay between images with 
the same winding number on opposite sides, we generally 
find a value which is much smaller than the time delay 
between images with different winding number. 

Moreover, for physically reasonable values of the coeffi- 
cients a, b, which are all of order one, the second term at 
the right hand side of Eq. (|41|l is much smaller than the 
first. For example, in the Schwarzschild black hole, the 
time delay between the first and the second relativistic 
images is (in Schwarzschild units) 



AT 2 



16.57 



(43) 



where the second term contributes only for 1.4% to the 
total time delay. 

For spherically symmetric metrics, we also have the 
very important relation 



(44) 



Namely, the dominant term in the time delay is not a 
new independent combination of the black hole metric 
function and gives no further hint for the classification 
of the black hole. On the contrary, the subdominant 
term is an independent combination and could be used 
in principle to constrain the black hole model. However, 
the subdominant term would be typically hidden below 
the observational precision and it becomes reasonable to 
approximate the time delay by its dominant contribution. 
In this way, a very interesting surprise arises. In fact, 
suppose we are able to measure the time delay between 
the first two images. Once we restore physical units, the 
ratio between this time delay and the minimum impact 
angle is 



AT 2 ,i „ D OL 
a = 2w ■ 



(45) 



where cq is the speed of light. In principle, by this for- 
mula, we can get a very accurate estimate for the distance 
of the black hole and hence of the whole hosting galaxy. 
The feasibility of such an estimate, will be discussed in 
the next section. 

In Appendix A, we treat the Kerr metric as an example 
of non-spherically symmetric black hole. In that case, 
most of the simplifications we have done, do not apply. 



IV. TIME DELAY IN SUPERMASSIVE BLACK 
HOLE LENSING 

In order to achieve a complete reconstruction of the 
characteristics of the black hole by strong field gravita- 
tional lensing, we must distinguish at least the outermost 



Local Group 


Mass 


Distance 


Schwarzschild 


Galaxy 


(M Q ) 


(Mpc) 


AT 2 ,i 


Milky Way 


2.8 x 10 6 


0.0085 


0.1 h 


NGC0221 (M32) 


3.4 x 10 6 


0.7 


0.2 h 


NGC0224 (M31) 


3.0 x 10 7 


0.7 


1.4 h 



TABLE I: Estimates for the time delay for the supermassive 
black hole located at the centers of three galaxies in the case 
of Schwarzschild spacetime geometry. The masses and the 
distances are taken from Richstone et al. 123. 



relativistic image from the others. Yet, as noted in jlflj . 
in order to achieve this, we need an optical resolution one 
or two orders of magnitude better than that reachable 
by short-term VLBI projects [22j. Therefore, relativistic 
images will possibly become a target for next generation 
projects. With this in mind, we can proceed to give es- 
timates for time delays between the first and the second 
relativistic images in realistic situations. 

We treat only black holes with spherical symmetry, be- 
cause only in this case we have the formation at first of 
more than one observable image. In fact, as noted in |17| . 
the phenomenology of spinning black holes is quite dif- 
ferent. In particular, if the source is not inside a caustic, 
only one image should become visible, while all the others 
stay very faint. On the contrary, if the source is inside 
a caustic, two additional non-equatorial images should 
appear. But an analytical treatment for these additional 
images is not available at present. 

Of course, we implicitly assume that the source must 
have temporal variations, otherwise there is no time de- 
lay to measure. Thus, an essential condition is that the 
source must be somehow variable. However, this is not a 
so restrictive requirement, since variable stars are gener- 
ally abundant in all galaxies. 

In Table m we present the values of the time delay for 
the black hole located at the center of the Milky Way and 
in other two galaxies of the Local Group. The results are 
obtained using the Schwarzschild metric. It is clear that 
we have a little chance to observe such short time delays 
for reasonable times of exposure. 

In order to have higher time delays, we need black 
holes with larger Schwarzschild radii, i.e. more massive 
black holes. At the same time we require that the mag- 
nification of the images must remain of the same order. 
Since we know that the magnification is proportional to 
A^Lens/-DoL, our request can be fulfilled if we consider 
lenses with a mass of two or three orders of magnitude 
larger than the black hole in the center of our galaxy, 
and located not farther than three orders of magnitude 
its distance. In this case, the measurement of the time 
delay becomes more favorable as shown in Tabled where 
we report our estimates for the time delay due to super- 
massive black holes located at the centers of not too far 
galaxies, according to spacetime geometry. 

The time delays range from few hours to several days. 
It must be kept in mind that a very deep exposure is 
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Mass 


Distance 


Schwarzschild 


Galaxy 


(M Q ) 


(Mpc) 


AT 2 ,i 


NGC4486 (M87) 


3.3 x 10 9 


15.3 


149.3 h 


NGC3115 


2.0 x 10 9 


8.4 


90.5 h 


NGC4374 (M84) 


1.4 x 10 9 


15.3 


63.3 h 


NGC4594 


1.0 x 10 9 


9.2 


45.2 h 


NGC4486B (M104) 


5.7 x 10 8 


15.3 


25.8 h 


NGC4261 


4.5 x 10 8 


27.4 


20.4 h 


NGC7052 


3.3 x 10 8 


58.7 


14.9 h 


NGC4342 (IC3256) 


3.0 x 10 8 


15.3 


13.6 h 


NGC3377 


1.8 x 10 s 


9.9 


8.1 h 



TABLE II: Estimates for the time delay for supermassive 
black holes located at the center of several nearby galaxies 
in the case of Schwarzschild spacetime geometry. The masses 
and the distances of the central black holes are taken from 
Richstone et al. |23| . 



needed to detect the very faint relativistic images. The 
precise time will depend on the characteristics of the fu- 
ture interferometers which will catch the relativistic im- 
ages and on the power of the source. However, we can 
imagine that an exposure of 10 hours can be still taken 
as a reasonable reference value for a deep imaging of a 
supermassive black hole. Then, with a high enough sam- 
pling and a suitable periodicity for the variable source, 
we can imagine to determine the time delay with an ac- 
curacy of few hours. So, most of the black holes in Tab. 
[n] would yield measurable time delays. 

Now consider the supermassive black hole in M87 and 
suppose we manage to reach an accuracy of 5% in a time 
delay measure. The resolution needed to resolve the first 
two images is 0.01 /zarcsecs, while the minimum angle is 
9 m = 11 //arcsecs. From formula (|45p. we can get the 
distance to M87 with an accuracy of 5% (the error in the 
angle measurement is negligible). This is already better 
than standard estimates by classical distance indicators 
|24| . whose accuracy ranges from 10% to 25%. 

So, gravitational lensing in the strong field limit may 
become a potentially competitive distance estimator in a 
not so far future. This is a consequence of the fact that 
the time delay is a dimensional variable and thus imme- 
diately leads to the measure of a scale. In the strong 
field frame, it is proportional to the mass of the black 
hole through the minimum impact parameter. However, 
it happens that we can also measure the minimum im- 
pact angle m — u m /DoL directly, so that their ratio 
leaves us with the distance to the lens. A time measure- 
ment can be done with a high accuracy and has the ad- 
vantage of being completely immune from any unwanted 
bias or systematics, unlike the classical estimates relying 
on luminosity measurements and typically highly model- 
dependent assumptions. 

A measurement as simple as this cannot be realized 
in weak field gravitational lensing, because it requires an 
accurate modeling of the gravitational potential. More- 



over, the length scale it measures is in general a more 
involved combination of all geometrical distances (-Dol, 
D ls , Dqs)- 

Another interesting possibility of strong field gravita- 
tional lensing is the possibility of getting a mass estimate. 
By the characteristics of the first two relativistic images, 
we can get the coefficients a and b, according to the pro- 
cedure described in [15j. They are generally sufficient to 
identify the class of the specific black hole. Afterwards, 
we can guess the theoretical u m in Schwarzschild radii 
for the specific black hole model. Combining with the 
the observed 6 m and with the Do l obtained by time de- 
lay, we get the Schwarzschild radius and hence the mass 
of the black hole. So, in principle, the time delay mea- 
surement would make the strong field gravitational lens- 
ing completely autonomous from external inputs coming 
from other methods. 

One final consideration about the subdominant term 
neglected in (|45|) : if we simply identify the time delay 
with its dominant contribution, we overestimate it by 1 
or 2%. However, once we have identified the black hole 
class by the coefficients a and b, we can easily evalu- 
ate the expected contribution of the subdominant term 
on the specific black hole model and subtract it from 
the observed time delay. We are then left with the pure 
dominant term and no more systematic errors (however 
small) are present. 



V. CONCLUSIONS 

Gravitational lensing in the strong field limit may rep- 
resent a key tool for the investigation of supermassive 
black holes. In principle, a complete characterization of 
the parameters of a black hole can be achieved by the 
study of the images formed by gravitational lensing of 
a background source. Technically, this study requires 
resolutions one or two orders of magnitudes better than 
actual VLBI projects, so that it stands as a possible ob- 
servational target for the next future. 

In this work we have pointed out that photons con- 
tributing to different strong field images take different 
times to reach the observer. This time delay is of or- 
der of few seconds for the black hole at the center of our 
Galaxy, but amounts to several days for more massive 
black holes at the centers of nearby galaxies. 

If the background source is characterized by an intrin- 
sic variability, it would then be possible to measure the 
time delay between different strong field images, with the 
important advantage of gaining a dimensional measure- 
ment for the scale of the system. This measurement can 
be immediately used to get an accurate distance deter- 
mination for the observed black hole, not affected by any 
kind of bias or model-dependent assumption. Identify- 
ing the black hole class by the use of the other strong 
field limit observables, we can also derive the mass of the 
black hole in a completely independent way. This result 
encourages our belief that gravitational lensing in the 
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strong field limit stands as an interesting (maybe pow- 
erful) method for the classification of black holes and 
the determination of their characteristics. Moreover, a 
new independent distance determination method is al- 
ways welcome in cosmological contexts. 
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hole, then we have to take care of the fact that one image 
will be direct and the other will be retrograde. We then 
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APPENDIX A: TIME DELAY IN KERR BLACK 
HOLES 

If the black hole is not spherically symmetric, the sim- 
plifications described at the end of Sect. IHII do not apply. 
In particular, the third integral in (|3 1|) does not identi- 
cally vanish. In this appendix we work out the time de- 
lay in a Kerr back hole as an example of non-spherically 
symmetric black metric. 

The Kerr metric projected on the equatorial plane 
reads 

A(x) = 1 - - (Al) 
x 

B ^ = i-Le (A2) 

X X 2 

C(x) =x 2 + a 2 + — (A3) 

x 

D(x) = 2~ (A4) 
x 

where a is the specific angular momentum of the black 
hole. 

We start directly from Eq. (|31[) . but we can still ex- 
press the dominant term AT in a simpler form. Con- 
sider first the case of two images on the same side of the 
black hole. Then everything works in the same way as 
for spherically symmetric black holes and we get 

A2£ ro = 27r(n-m) = . (A5) 
a 

Of course, the values of a and a depend on the sign of 
the spin a, i.e. they are different for photons winding in 
the same sense of the black hole (direct photons) and for 
photons winding in the opposite sense (retrograde pho- 
tons). 

If we wish to evaluate the time delay between two rel- 
ativistic images appearing on opposite sides of the black 



FIG. 2: Time delay (as a function of the black hole spin) 
between the second and the first relativistic images appearing 
on the same side of a Kerr black hole. 

get 

AT° m = ^ [2™ + 7 - b(a)] + b(a) - 
a(a) 

|M[27rm-7-6(-o)]-6(-o) (A6) 

and we see that now we also need the coefficients b and 
b for the calculation, since they are not the same for 
the two images and do not cancel like in the spherically 
symmetric case. 

Considering a source aligned behind the black hole 
(7 = 0), in Fig. |2]we plot the time delay between the 
second and first images appearing on the same side of 
the black hole. For positive a the two images are direct 
and for negative a they are both retrograde. We see that 
the time delay decreases if the images are both direct, 
while increases if they are both retrograde. We can also 
notice that the largest contribution to the time delay still 
comes from AT, while the second term in l|31|) at most 
contributes for 6% when a = 0.5 and the last term stays 
below 0.7%. 

The situation is quite different for images on opposite 
sides (Fig. The time delay is zero when a = and 
becomes negative for positive a. This means that direct 
light rays take less time than retrograde rays to wind 
around the black hole. This is naturally understood since 
the radius of the photon sphere is larger for retrograde 
light rays. For high values of the black hole spin, AT-j^ 
becomes comparable to ATf-^ Another interesting fact 
is that the second term in (|31|) is of the same order of the 
dominant term AT, the ratio being roughly —1/3. The 
last term also contributes for 1.7%. 
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